Coulomb "blockade" of Nuclear Spin Relaxation in Quantum Dots 
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Wc study the mechanism of nuclear spin relaxation in quantum dots due to the electron exchange 
with 2D gas. We show that the nuclear spin relaxation rate T^ 1 is dramatically affected by the 
Coulomb blockade (CB) and can be controlled by gate voltage. In the case of strong spin-orbit (SO) 
coupling the relaxation rate is maximal in the CB valleys whereas for the weak SO coupling the 
maximum of 1/Ti is near the CB peaks. 
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Discreteness of energy spectra in quantum dots (QD), 
in combination with the Coulomb electron-electron in- 
teractions, determines all their physical properties [Q. 
In this paper, we study the impact of these two crucial 
features of QD on the interaction of electrons with nu- 
clei. In particular, we evaluate nuclear spin relaxation 
(NSR) rate, 1/Ti, which determines the effectiveness of 
dynamical nuclear spin polarization 0] as well. 

What mechanisms of NSR are efficient in the QD? Be- 
cause even weak magnetic fields suppress the NSR due to 
the dipole-dipole interaction of nuclear spins [^|, we will 
consider NSR caused by the hyperfine coupling (HC) be- 
tween nuclear, S, and electron, cr spins: 

H h = AVj2 s An), B{n) = 9+(T i )tr aP 9p(T i ), (1) 
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where A is the HC strength, and V is the crystal cell vol- 
ume. Substitution of the operator B by its mean value, B, 
leads to the Hamiltonian of the system of nuclear spins, 
St in an effective magnetic field, B. Being uniform, this 
field causes only a uniform precession of Sj, rather than 
NSR. Thermal fluctuations of B in metals lead to the Ko- 
rringa NSR mechanism || - each nuclear spin flip is ac- 
companied by the creation of a triplet electron-hole pair. 

In QD this process violates the energy conservation: 
the electron spectrum is discrete and Kramers degener- 
acy is lifted by the Overhauser nuclear field and/or ex- 
ternal magnetic field. In order to overcome the mismatch 
of electron and nuclear spin splittings, a level broadening 
Tq was included for electrons localized on donors 0,[|]6j 
and in the quantum Hall effect regime f]J|||. However, 
To must be introduced with a great caution. Indeed, any 
broadening implies that the final state of the system is 
not a discrete eigenstate of the QD but rather an eigen- 
state of a larger system with a continuous spectrum which 
includes QD as its subsystem. Thus, one has to specify 
the nature of the QD coupling to the outside world. 

We propose that charge exchange between the QD 
and the reservoirs (with continuous spectrum) leads to 



NSR. Our main result is that NSR can be controlled and 
changed by orders of magnitude by the gate voltage, V g , 
which defines the average number of the electrons in the 
QD. The origin of this effect is the Coulomb blockade 
(CB) [0 that determines the tunneling of electrons in 
and out of the QD. Furthermore, l/Tj in QD is governed 
not only by the QD electrons tunneling rates and the ef- 
fective fields B(vi), but also by the probabilities for the 
QD to have a particular integer charge, Q, as well as by 
the occupation numbers of the electronic states in the 
QD and in the leads at a given Q. Moreover, we find 
that the NSR rate depends on the sensitivity of the QD 
spectrum to the change in the nuclear spin configuration, 
which is maximal in the presence of electron spin-orbit 
coupling (SOC). 

Let us now discuss the role that SOC plays in NSR. 
Due to SOC the total spin of nuclei and electrons is not 
conserved. Each electron level, even if Kramers degener- 
acy is lifted, contains mixture of up and down spin states, 
and energy conservation law does not prevent nuclear 
spin flips (nuclear spin splitting Hlu — > 0). Does the pres- 
ence of SOC leads to NSR even in absolutely closed dots? 
The answer to this question is negative for the following 
reason. Due to SOC, both the magnitude, and the direc- 
tion of the effective magnetic field, B =< B >, become 
spatially inhomogenous but remain time-independent. 

Although the field originating from SOC, B — £> so (r), 
causes an inhomogeneous spin precession and thus broad- 
ens the NMR line, it can not lead to the actual NSR. The 
reason is that B so (r) is time- independent. One can sepa- 
rate this inhomogeneous spin precession by the spin echo 
technique [Q. Furthermore, a time-independent B can re- 
duce the nuclear polarization but is unable to eliminate 
it altogether. For example, if local nuclear polarization 
is parallel to B, it will not change at all. 

However, even weak electron tunneling on and off the 
dot causes true NSR. The reason is that B so is determined 
by the particular electronic configuration, and is modified 
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by each tunneling event, thus becoming time- dependent. 
This dynamics leads to NSR, because it is not reversible. 
Indeed, due to inhomogeneous precession of the nuclear 
spins, the adiabatic electron eigenstates change between 
tunneling events, and it is improbable that the system 
returns to the same electron/nuclear configuration. 

Let us now roughly estimate the NSR rate due to 
this mechanism, given the effective electron escape rate 
7. Typical rotation angle of S within the time 7 _1 is 
B so /j -C 1. At larger time scale, t, nuclear spins undergo 
random spin diffusion. The mean square angle deviation 
is (Q 2 (t)) ~ (B so /T) 2 tj = B 2 so t/~f. Defining NSR rate 
1/Ti by (6 2 (Ti)) ~ I, we obtain 1/Ti - B 2 Jj. The 
rate 7 can be fine tuned by adjusting the gate voltage, 
V g , applied to the QD. 7 reaches maximum at the CB 
peak, where the dot conductance G is maximal, and is 
minimal in the CB valley. We reach a counterintuitive 
conclusion: NSR rate peaks in the CB valleys, while near 
the CB peaks it is suppressed, see Fig. 2. As we shall see, 
this conclusion holds as long as SOC is not too weak. 

The described physical situation is not unique for 
NSR in QD. In particular, it resembles Mandelstam- 
Leontovich-Pollak-Geballe relaxation mechanisms of ab- 
sorption of infrared radiation jL2|: As long as the relax- 
ation rate 7 r exceeds the frequency of the radiation to, the 
relaxation rate is inverse proportional to 7,. . One can also 
recall the motional narrowing of spectral lines fjlj ]. 

Our goal is the microscopic theory of NSR in QD 
which includes the calculation of 7 and B so . We start 
with the standard Hamiltonian of a QD connected to 
leads (see e.g. Ref. |l| for the discussion of its valid- 
ity) H = H D + H L + H LD , where H D (t), H L and H LD 
describe, respectively, the dot, the leads and tunneling 
between them: 



H 



D 



(2) 



HL = T,k,j e k,j b k,j b hj>H LD = Ek,j;a*£ 6 fe,4 ffl a + h - C " 

Here a^a^ and bk,j,,b' k . are creation/ annihilation elec- 
tron operators describing single-electron QD states a, 
and electrons in state k in the lead j, respectively; 
h = J2 a a a a a is the operator of the number of particles, 
N = CV g /e, C is the capacitance of the dot, e is the elec- 
tron charge, and E c = e 2 /2C is the single-electron charg- 
ing energy; matrix elements P a determine the tunneling 
widths of QD states \a), 1$ = J^j ^{t^u, where v is 
the density of states in leads. 

The NSR rate 1/Ti due to HC is connected to the lin- 
ear response of the electronic system to Zeeman magnetic 
field induced by nuclei at their location M : 



Tf 1 = 2A 2 V 2 TIm Xxx (u)/u; 



(3) 



The transverse spin susceptibility \xx is defined by the 
response 8{s x (r)) u = Xxx(u)hu, h u = J dtexp (iujt)h(t), 
of the electron spin density, s x (r), to the perturbation 



6H = h(t) V [sxirJlapalap. 

* — ' an 



(4) 



As soon as SOC is included, the diagonal matrix elements 
[sa;(r)] QQ , become finite and result in the shifts of the en- 
ergy levels e a . We assume the dot to be weakly coupled 
to leads which are in thermodynamic equilibrium. Then 
the probability, Pj, for the QD electrons to be in the 
state j, tends to follow the Gibbs distribution Pj 3 deter- 
mined by the instantaneous value of H(t). However, Pj 
is delayed as compared with H(t) due to finite relaxation: 



P 3 {t) = Pf{H{t)} 



(5) 
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The matrix T characterizes the relaxation of the devia- 
tion of Pj from their equilibrium values: 



dt 



Pi(t) = -Y,TihPk, £ r '-* = E I v , * G = - (6) 



The second term in the RHS of Eq. gives rise to the 
relaxational dissipation |L2] . The average spin equals to 



Sx(r)) = ^ n J a PjS a , s a = [s x (r)] aa 

* — 'ja 
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where nP a — 0, 1 is the occupation number of one-electron 
level a for a given state of the electron system j. 



a) 




b) 




FIG. 1. Processes dominating NSR in the (a) vicinity of 
and (b) away from the Coulomb blockade peaks. Lines de- 
note single-electron levels. Dots indicate filled levels. Arrows 
show the relevant transition. Dashed lines in (b) show the 
actual electron processes leading to the relevant transition. 

Substituting P 3 (t) into Eqs. (0) and Eq. (g), 

we obtain 

Xxx{u) and Tf 1 : 



1 = _ A 2y2 T y 

Ti ^. 



dhe'" 1 



dPf 



We emphasize that this NSR rate is inversely propor- 
tional to the rate of the population relaxation of the QD 
states, provided the latter rate is larger than u>. 

As follows from Eq. (JsJ) , the NSR is dominated by elec- 
tron configurations such that (i) their equilibrium prob- 
abilities are most sensitive to the shifts of one-electron 
levels, and (ii) their relaxation is slowest. Let us iden- 
tify the optimal configurations for temperatures, T, be- 
low the mean level spacing, A. We begin with V g 
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tuned to a CB peak. The energy distance to the peak 
[/(TV) = E C {N ~ N ) ~ A, with N being a half-integer 
closest to Af. Then the relaxation occurs (see Fig. [l]a) via 
the two-level system (TLS) of the QD states with electron 
occupation numbers A^ ± 1/2, which we denote 1 and 
respectively. Occupation of other states is suppressed ex- 
ponentially, i.e. Pq + Pi = 1. The Gibbs distribution for 
such TLS gives P G /P G = exp [-(U(Af) + ei)/T] where 
€i is the energy of the filled electron level in state "1". 
The relaxation matrix for the two states is determined 
by single tunneling width Y. 



i) 
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where / = {1 + exp [(U(Af) + ei)/T]} _1 is the electron 
occupation number in the leads. From Eq. (ph, we find 
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(10) 



Thus, NSR decreases rapidly as V g deviates from the CB 
peaks, because the sensitivity of Pj to e a is exponentially 
suppressed by lifting the charge degeneracy. 

Naively, one may expect the Eq. (|l0|) to hold, at least 
qualitatively, even far from the peaks. However, as soon 
as \U(Af)\ exceeds A, another TLS becomes optimal, see 
Fig. [Lb. The energy difference between the two config- 
urations in this TLS (we will call them 1 and 2) is e e h 
- the smallest energy of the electron-hole pair excited in 
the dot. At \U(Af)\ > A, e eh < \U(Af)\, the population 
of all other states is negligible. As a result, Pi + P 2 = 1. 
In the equilibrium, P 2 G /P G = exp [-(e 2 (t) - ei(f))/T], 
i.e. in contrast with the previous case Pi )2 are \U(N)\- 
independent because now the relevant excitation does not 
change the charge of the QD. Charging energy, however, 
determines the rate of the transitions between the two 
states. If V g is not too far from the peak, U(Af) < 
Tln(.E^/IVA), the transition 1 — > 2 proceeds via the 
state (i.e. due to electron transfer between the QD and 
the reservoir) as an intermediate real state. One finds 
EOT 
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(11a) 



This result can be derived from the master equation in- 
cluding the state 0. For its interpretation consider the 
rate of the 1 — > 2 transition, [T]n- The electron es- 

cape rate out of the state 1 is 1^ e t . After the 
electron leaves the QD, another one enters during the 
time ~ 1 /r t . It occupies the state 2 with the probability 



The product of those two factors gives 



p(2) /(p(l) , t-(2) 

the probability of the transition 1 — ► 2. 

At larger deviations from the CB peak, |{7(Af)| > 
Tln(£'^/r t A), processes changing the QD charge can be 
neglected. The relaxation is now determined by the in- 
elastic co-tunneling mechanism JTl|,p| , i.e., the state is 
used as a virtual state. The Eq. (11a) becomes 
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where £21 = £2 — £1 > 0, is the lowest energy of the 
electron-hole excitation in the QD, F = coth(e2i/2T). 
The factors in parenthesis arc determined by the phase 
volume for the electron-hole pair created in the reservoir 
as the result of the inelastic co-tunneling process. 

Using the equilibrium occupation numbers of the two 
QD states and Eqs. (lib) and (0), we find 



r t - 



A 2 V 2 \{s(r) x }ii - {s(v) x } 22 \< 
Ah max (r t ,r c )cosh 2 Si 



(12) 
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This result qualitatively differs from Eq. (|10|), and the 
dependence of 1/Ti on V g is non-trivial, see Fig. ||. 
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FIG. 2. Sketch of the predicted NSR dependence on V g in 
the presence of SOC. Dashed line shows the linear conduc- 
tance through the dot. 

With yet further change of V g towards the bottom of a 
CB valley, one has to take into account the state with 
an electron removed from the dot (hole-like process) and 
with an electron added to the dot (electron- like process) . 
This changes the rates r tjC in Eq. (|12f): 



4r^ 1) rp ) cosh ^e~ Ec / T \U(M)\ 
It = — ^ cosh 
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(13) 



Note that the inelastic rate T c vanishes at the bottom of 
the valley, because the electron- and hole-like processes 
with the same final states are coherent, and their am- 
plitudes are opposite in sign due to the Fermi statistics 
(the higher levels would give non-zero contribution, but 
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the conclusion about the maximum of 1/Ti at the bottom 
of the CB valley remains valid) . Overall dependence of 
1/Ti on V g in the presence of SOC is sketched on Fig. ||. 

We now discuss NSR when SOC of electrons in QD 
is absent, and each QD level corresponds to a spin pro- 
jection aligned or antialigned with magnetic field. For 
typical distance between the lowest unoccupied QD level 
level and highest occupied QD level with opposite spin, 
£2 - £i ~ A > , T, the calculation of the suscepti- 
bility Xxx(w) is equivalent to the evaluation of the Fermi 
golden rule probability of the simultaneous electron and 
nuclear spin flip. In this process, the initial and final 
electron states are real, one or both of them are in the 
reservoir, and energy-nonconserving electron spin flip in- 
side the QD is incorporated as the virtual transition. 

In the vicinity of the CB peaks, the two relevant QD 
states are electron on the lower level (1) and electron in 
the leads (2), see Fig. 0. Deeper in the valleys the main 
contribution is due to the process analogous to the elastic 
cotunneling in the CB valley conductance (Fig. ||b): the 
initial and final QD electron configurations coincide but 
the electron-hole triplet pair with is excited in the leads. 



a) 




FIG. 3. Dominant NSR processes at weak SOC in the (a) 
vicinity of and (b) away from the CB peaks (e% - ti > T^). 
Solid lines with arrows denote transitions between initial and 
final states, dashed lines with arrows correspond to stages of 
those transitions. The spin flips occur in transitions "A". 

The result for the overall NSR rate dependence on V g in 
the absence of SOC is 
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Such NSR rate [14) is directly proportional to the rate of 
the population relaxation of the QD states, is maximal 
at CB peaks, minimal in the valley, and much smaller 
than 1/Ti in the presence of SOC. 

We now estimate the NSR rate. In GaAs, A — 40/ieV 
H , and V — 45 A 3 . We consider QD formed in a 50A-wide 
quantum well, with 0.1/Ltm in-plane dimensions, charging 
energy E c — 0.5meV, level spacing A = 50^eV, and tun- 
neling width r t = 2/ieV, at T = lOOmK. The typical 
matrix element sn is of the order of the inverse volume 
of the dot. Therefore we use V 2 s\ 1 = 0.8 x 10 -12 in 
, and obtain T, 
2 x 10 

CB peak, at U{N) = 55/xeV, Fig.2. In the CB valley, by 
using Eq. (n3il3), and assuming su — S22 ~ O.lsn, at 



Eq. (|1C 

Tf 1 



0.8 x 10" 12 
! — 0.25Hz at the CB peak, and 
3 Hz in the NSR minimum close to the 



[/(W) = 0.495meV we obtain Tf 1 = 1.6 x 10 3 Hz. Our 
analysis shows also that even for QD with rather weak 
SOC, Tf 1 in the CB valley has at least local maximum, 
so that the minimal NSR rate occurs at V g values be- 
tween the CB peak and the CB valley bottom. Thus, V g 
indeed changes Tj~ by orders of magnitude. 

In conclusion, we have shown that nuclear spin relax- 
ation in QD is strongly affected by the Coulomb block- 
ade, as well as the spin-orbit coupling, and proposed the 
relaxational mechanism of NSR. The NSR rate is pre- 
dicted to have non-trivial dependence on V g . Similarly 
to nuclei in experiments on quantum point contacts |l5[ , 
nuclei shall affect transport properties of QD. The gate 
voltage depedence of NSR can be used, in particular, for 
creating and sustaining nuclear spin polarization in QD 
for electron spin filtering (we are grateful to CM. Mar- 
cus for discussions of this experimental setting). This 
work was supported by DARPA QUIST program (Y.L.- 
G.) and (B.A.) and by Packard Foundation (LA.). 
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